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Abstract 

The exact equation of spin motion in a cylindrical coordinate system with allowance for electric 
dipole moments of particles has been derived. This equation is convenient for analytical calculations 
of spin dynamics in circular storage rings when the configuration of main fields is simple enough. 
The generalized formula for the influence of a vertical betatron oscillation on the angular velocity of 
spin rotation has been found. This formula agrees with the previously obtained result and contains 
an additional oscillatory term that can be used for fitting. The relative importance of terms in the 
equation of spin motion is discussed. 

PACS numbers: 12.20.Ds, 29.27.Hj, 41.75.-i 
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I. INTRODUCTION 



New experiments for measurement of the muon anomalous magnetic moment (AMM) 
and the electric dipole moments (EDMs) of fundamental particles (3, \$ in storage 
rings need an extremely high accuracy of measurement of the spin precession in storage 
rings. To distinguish the spin precession induced by the EDMs of particles and other minor 
effects, the spin motion has to be mathematically described in a convenient manner. The 
goal of this paper is to derive the equation in cylindrical coordinates which allows an exact 
analytical description of spin dynamics in storage rings. We do not confine ourselves to 
specific configurations of electric and magnetic fields. In particular, electric or magnetic 
focusing fields and an accelerating electric field can be used. The influence of EDMs on 
the spin dynamics is taken into account. We consider the problem of averaging the angular 
velocity of spin rotation and calculate the generalized formula for the influence of vertical 
betatron oscillation upon this quantity. 

There are many algorithms of computer calculations based on different coordinate sys- 
tems (e.g., on Frenet-Serret coordinates ||). These algorithms cover any problem of beam 



and spin dynamics in storage rings. However, an analytical solution of the problem can 
)e necessary for several high-precision experiments. We mean the g— 2 Q, 0] and EDM 
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experiments, those precision is extremely high. These experiments need a clear un- 
derstanding of spin dynamics. The use of cylindrical coordinates for analytical calculations 
of spin dynamics can be very successful if the configuration of main fields is simple enough. 
For circular storage rings, it is quite natural to describe the spin motion just in cylindrical 
coordinates. For this reason, other formalisms are less convenient for the g— 2 and EDM 
experiments. This problem will be considered properly later. 
The relativistic system of units h = c = 1 is used. 



II. GENERAL EQUATIONS OF PARTICLE AND SPIN MOTION 

The particle trajectory in a circular storage ring is approximately a circle. Coherent 
betatron oscillations (CBOs) of a beam as a whole in the horizontal and vertical planes affect 
the trajectories of individual particles. Synchrotron motion occurs whenever radio frequency 
(rf) cavities are used. The incoherent motion of particles can also take place. As a result, the 



2 



particle trajectories are never closed. Field defects and misalignments of magnets, electric 
field plates and electrostatic quadrupoles cause particle trajectory distortions. These effects 
lead to a change in the magnetic field acting on the spin in the particle rest frame. As a result, 
the spin motion becomes complicated. For this reason, equations of spin motion correctly 
taking into account perturbations of the particle trajectory have only been obtained for some 
specific problems jslQ, In Refs. ^L,\^, the corrections to the spin rotation frequency due 
to the vertical and radial betatron oscillations have been calculated. 

As a rule, the use of the one-particle approximation is sufficient to describe the particle 
and spin motion. In this approximation, coherent and incoherent betatron oscillations give 
similar effects. However, the CBOs are more important than the incoherent BOs because 
they cause a systematic shift in the spin precession frequency in the real experiments (see 
Refs. 

The particle motion is described by the Lorentz equation 

^ = e(E + f3xB), (3=^=P-. (1) 

dt c 7m 

It is convenient to use the unit vector of momentum direction N = p/p which defines 
the direction of particle motion. Since 

dN p p 

dt p p 3 

Eq. (JTJ) takes the form 

dN _ e / N x E\ 

wxJV, uj = I B — , (2) 



dt 7m \ (3 

where uj is the angular velocity of particle rotation. 

The spin motion is determined by the Thomas-Bargmann-Michel-Telegdi (T-BMT) equa- 
tion 

ds 

— — Ht-BMT x s, 

dt 

n rf 

where s is the spin vector. This equation was derived by Thomas l£ ] (also by Frenkel [1CJ) 
and, in a more general form, by Bargmann, Michel and Telegdi 
means an expectation value of a quantum mechanical spin operator. 



. In this work, "spin" 
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The T-BMT equation describes the motion of spin in the rest frame of the particle, 
wherein the spin is described by a three-component vector. The longitudinal direction in 
this frame is mapped from the longitudinal direction in the laboratory frame. 

A comparison of Eqs. (J2J) and (j3J) shows that the spin of the particle placed in a vertical 
magnetic field rotates in the horizontal plane with a frequency proportional to g— 2. We use 
the term "g— 2 precession" for any spin rotation in the horizontal plane, even in the presence 
of an electric field. 

III. CORRECTIONS TO THE ANGULAR VELOCITY OF PARTICLE MOTION 
IN THE HORIZONTAL PLANE 

The particle spin motion in storage rings is usually specified with respect to the particle 
trajectory. The use of the cylindrical coordinates considerably simplifies the analysis of g— 2 
precession and effects caused by the particle EDM. The particle rest frame is less convenient 
for taking into account the acceleration of the particle and the tilt of its orbit. When 
the configuration of the main fields is simple enough, the use of cylindrical coordinates 
simplifies the form of equations of particle and spin motion. However, the description of 
the spin motion in the cylindrical coordinate system is a difficult problem because the axes 
of this system are defined by the position of the particle being in oscillatory motion. The 
transformation of the T-BMT equation to the cylindrical coordinates should be performed 
with allowance for the oscillatory terms in the particle motion equation. 

The vertical CBO and orbit distortions change the plane of particle motion. The 
(pseudo)vector of angular velocity becomes tilted. The instantaneous plane of particle mo- 
tion does not coincide with the horizontal plane, and the instantaneous plane of rotation 
of vector N is not horizontal. The angle $ between two positions of the rotating vector 
N in the tilted plane is not equal to the angle between two corresponding horizontal 
projections. Therefore, the vertical CBO and the orbit distortions change the instantaneous 
angular velocity of particle motion. This effect can be calculated. 

We suppose the x- and y-axes are horizontal and the z-axis is vertical. Since we use 
the cylindrical coordinate system, it is convenient to direct the z-axis orthogonally to the 
plane of the unperturbed particle motion. We can define the angle of particle rotation in 
the x?/-plane as the angle between two horizontal projections of the vector of momentum 
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direction in two indicated positions, N\\ and N». The infinitesimal angle of particle rotation 
in the xy-plane, dtp, is given by 

(JVii x M) ■ e z fJVn x dNn) ■ e, 
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|iV|||-|JVj| |JV||, 

where dN\\ = N» — JVy and dJVy is an infinitesimal vector. The sign || means a horizon- 
tal projection for any vector. The quantity dNn determines the deflection of the particle 
momentum direction in a time interval dt. The instantaneous angular velocity of particle 
rotation in the horizontal plane is equal to 

• _ dtp _ (JV|| x N\\) ■ e z 
(/)= H = \N] 

where 
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u z - o, (4) 



{u x N x + LO y Ny)N z _ {u p N p + u)fNt)N ± 



1-iV 2 i-^ 2 2 

Different components of vector uj are defined by Eq. (J2J). Indexes p and mean projections 
onto the basis vectors e p and of the cylindrical coordinate system. 

Eqs. (JU),© ar e exact. We can confirm the validity of these equations by the fact that 
they perfectly describe a constant tilt of the particle orbit and then evaluate the quantity o 
in experimental conditions when the tilt of the particle orbit oscillates. 

If the particle orbit is perfectly horizontal, the momentum direction is defined by 

N x = -sm(u t + (p ), N y = cos (w t + ipo), N z = 0, 

where uq is the cyclotron frequency and <^o is an arbitrary phase. If the normal to the tilted 
particle orbit is orthogonal to the y-axis and deflected from the z-axis by a constant angle 
9, y-components of all vectors remain unchanged and the components of vectors uj and N 
are equal to 

ui x = ujq sin 9, ujy = 0, uj z = u cos 9, 

(6) 

N x = — cos 9 sin (u t + (p ), N y = cos (uot + ipo), N z = sin 9 sin (u t + (p ). 



In this case, Eq. (@J takes the form 

• u cos 9 



(7) 



1 — sin 2 9 sin 2 (uj t + <p ) 
As a result of integrating and averaging over time, we obtain < <fi >= uq. This result 
agrees with the evident fact that the average frequencies of particle motion in the tilted and 
horizontal planes are equal and therefore confirms the validity of Eqs. P |L (j5 |) . 



To evaluate the quantity o in real experimental conditions, we may confine ourselves 
to taking into account only the particle trajectory perturbations caused by the radial and 
vertical CBOs. The synchrotron motion does not affect this quantity. By our estimate, it is 
sufficient to approximate the horizontal and vertical betatron motion by the simple forms 

N p = — = p sin (uj r t + a), 
V 

N z = — = i>Q sin (uj v t + 5), 
V 

where po an d V'o are angular amplitudes, u r and u v are angular frequencies of the radial and 
vertical CBOs, respectively. 

With allowance for orders of quantities 

u p ^ip u v , 10$ ~ palpal, iV p ~ p , « ±1, N z ~ ip , (8) 

we obtain that the quantity o is of the third order in the angular amplitudes po an d tpo- 
Moreover, it oscillates and therefore averages to zero. If we take into account only the 
second-order terms in the angular amplitudes and the average particle orbit is not tilted, 
the quantity o is negligible. Approximately, 

* = u . = -±( B .-(»2L*k). (9) 

IV. EQUATION OF SPIN MOTION IN CYLINDRICAL COORDINATES 

To transform the general equation of spin motion (jHJ) to cylindrical coordinates, we need 
to find the quantities 

ds p dsA) ds z 
——, — — , and — — . 
at at at 

It follows from the geometry of the problem that the horizontal axes e p and rotate with 
the instantaneous angular velocity 

u>' = 4>e z . 

It is easy to show by simple transformations that the spin motion with respect to the 
axes of the cylindrical coordinate system can be written in the form 

ds 

— = uj a x s, uj a = fL- 4>e z . (10) 

In this equation, uj a is the the angular velocity of spin rotation in cylindrical coordinates. 
The corresponding angular velocity in Cartesian coordinates equals f2. If the spin motion 



is described by the T-BMT equation, f2 = SIt-bmt- The difference between the quantities 
uj a and f2 is caused by the rotation of the axes e p and e^. 

The equation of spin motion with allowance for the EDM can be obtained by modifying 
the T-BMT equation. The EDM and AMM determine the real and imaginary parts of the 
same interaction Lagrangian 3, Q|- This Lagrangian is used to derive an extended Dirac 
equation. The contributions of the EDM and AMM to the Lagrangian are equal to 
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Camm = j^F^, C EDM = -i^a^F^, = % - (-f-f - 7 V) , (H) 

where // = e(g - 2) /(4m) is the AMM, d is the EDM, 7^ (/1 = 0,1,2,3) are the Dirac 
matrices, and F^ y is the electromagnetic field tensor. The use of the explicit forms of Camm 
and Cedm shows that the contribution of the EDM to the full Lagrangian can be obtained 
from the corresponding contribution of the AMM by the substitution 

B -> E, —B, p! -> d. (12) 

Takin g int o account the particle EDM leads to the following modification of the T-BMT 



equation 



ds 

— — ft X S, Q = Qt-BMT + ^BDMi 

dt 



Vedm = (e - ^7 W ■E)+f3xB ) j, (13) 

where £It-bmt is defined by Eq. (j2J) and 77 = 4dm/ e. 

The EDM affects the particle motion if the electric field is not uniform. However, the 

correction to the equation of particle motion is negligible. If rf cavities are not used, we can 

also neglect the term j3(/3 ■ E). 

7 + 1 

The transformation of Eq. (|13|) to cylindrical coordinates is given by Eq. (|lUp. where 

u a = -- Lb - ^-W • b) 

ml 7 + I 



1 \ 1 

a )(j3 x E) + 



7 1 / 7 



+ |(^-^ T W-^)+/3x^}+oe z , a= 9 -^. (14) 

This formula is exact, and u> a is the angular velocity of spin precession. Eqs. (JTUJ) , (JHJ) 
describe the spin motion in arbitrary storage rings with allowance for the EDM. The term 
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oe z in formula (|14|) is negligible for the EDM and g— 2 experiments. Measurements of the 
AMM in the g— 2 experiment are performed at an energy such that l/(7 2 — 1) = a, that is 
7 = 29.3. In this case, the third term in Eq. (|14|) is equal to zero Q, 2|. 

After neglecting the small terms, the equation for the angular velocity of the g— 2 pre- 
cession with allowance for the EDM takes the form 




/3(/9 -B) 



-a)(/3x E) 



1 

+- 

7 



7/ 



+ l(E + (3xB) 



(15) 



Formulae l|14 J) . (|15J) are useful for analytical calculations of spin dynamics in cylindrical 
coordinates with allowance for field misalignments and beam oscillations. 



V. AVERAGING OF THE ANGULAR VELOCITY OF SPIN ROTATION 



The spin dynamics given by formula (jH| is simple when the tilt of the spin out of the 
xy-plane accumulates from turn to turn. Such a behavior of the spin takes place in the 
planned EDM experiment. However, in the g— 2 experiment one needs to determine the 
average and integral characteristics of the spin motion in the horizontal plane 

HQ. 

The 

instantaneous angular velocity of spin rotation in the horizontal plane, ip, is characterized 
by the change of angle ip determining the spin orientation in this plane. The quantity ip can 
be found similarly to the corresponding quantity <fi describing the particle rotation (Section 
III). It is given by 



where 



O 



v dt 



M x £x + Myiy 



[OJ n 



O. 



6 
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(16) 



(17) 



and £ = s/s. 

As an example, we can derive the generalized formula for the average value of ip affected 
by the vertical CBO. When the particle and spin motion is unperturbed, the angular velocity 
of spin rotation is equal to 



\ea\ 



-B z . 



m 



(18) 



We suppose that B z > 0. 

The vertical CBO gives an important correction to the g— 2 frequency (the pitch correc- 
tion). The formula used currently for this correction has been calculated by Farley M|. T he 
result has been confirmed by Field and Fiorentini pj and by computer simulations [l7[. If 
we ignore the possible existence of the EDM, then Eq. (|15|) can be written in the form 

^ = {a + a 3 cos [2(u v t + 5)]} (e z x £) 
+a 2 sin (uj v t + 5) (e^ x £) + ai cos (a> v t + 5) (e p x £) , 

where 

a = Afi U - "^~^o 

7 — 1 

«i = A/c^o, «2 = -^o V'o, 

-1 7 

a 3 = Afi — — ^o> 

and A is equal to 1 and —1 for negative and positive muons, respectively. The factor / is 
given by 

/ = l + a7 _l±^ = l + a ^ 7 _I 

7 7 

and 

/ = 1 + a 7 

for electric and magnetic focusing, respectively. In the muon g— 2 experiment, ipo ~ 1CT 3 . It 
can be shown that the effect of 03 on the spin dynamics described by Eq. (jl9)l is negligible. 
If we substitute the zero-order approximation 

^ O) =£||cos(a o t + O ), 4°)=£||sin(a o t + o ), £<°> = £ ± (20) 
into the right side of Eq. (|19|) , the first-order approximation is given by 
£ p = f|| cos (aot + O ) + aoa ^ + ^v® 2 £ ± cos (u;„t + 5), 



05 - ^,1 



4ll sin (a t + 0o ) H 5 5— £± sin (u;„£ + 0), 



& = £j 



a 5 -< 

a a 2 + u v ai . . . 
^ — sm (a t + 0o ) sm (a;„t + dj 



(21) 



05 -w w 



a ai + ^a 2 , , , , . « 

H o 5 — cos (a t + 0o ) cos (a;„t + <)) 



where £j_ =< £ z (t) > and £|| = y 1 — Angular brackets mean the time average. If |a | is 
close to u v , the corrections to the spin dynamics given by Eq. (|2*T|) can be large. Therefore, 



the condition |eto| = oj v defines a resonance. When the spin rotation frequency is far from the 
resonance, Eqs. (fTBj) . (|T7|) . (fT^j) . (PT|) lead to the following expression for the average angular 
velocity of spin rotation: 



<ip>=\n (i-c F ), c F = l - 



% ^(/-l)(/-l + 2/ 7 ) 



K (22) 



This formula coincides with the result obtained in Refs. (3, 0| • However, in the real g— 2 
experiment u v /Qq ~ 10. Therefore, we may not average over the spin rotation which is 
rather slow. As a result of averaging only over the vertical betatron oscillation, Eq. (|2*2*J) 
takes the form 

< j, >= AOo(l - C), C = C F + (7 ~ ^"jj " cos P( n o* + A0oM " ^4- 

(23) 

Since the average value of the last term oscillating with the angular frequency 2Q is zero, 
formula (J2H|) is in the best agreement with the previously obtained result. However, it is 



reasonable to include the oscillatory term in the fitting process instead of eliminating it. In 
the real j 
given by 



7T 

the real g— 2 experiment, / ~ 1, (1 = 0, 7 > 1, 0o equals ± — , and the coefficient C is 



C= |[l-cos {2Q t)]^l 

The quantities Cp and C are rather small because they are proportional to i^q. In the real 
g— 2 experiment, ipo is less than 1 • 10~ 3 . 

The vertical betatron oscillation violates the sinusoidality of spin motion. 

Essentially, the result of averaging over the vertical betatron oscillation is independent 
of an initial phase of this oscillation. Therefore, coherent and incoherent oscillations give 
the same effect. The oscillation of the angular velocity < ip > averaged in such a way is the 
result of the fixed initial spin direction. 

The quantity O in Eqs. (jlfij) and (fTTjl is much smaller than the corresponding quantity 
o in Eqs. (J1J) and (jHJ). A significant difference between these quantities is conditioned by a 
relative smallness of N p and in comparison with £ p and (w a )^. Such a smallness is due 
to the fact that the radial component of particle momentum is proportional to the small 
quantity po because it is caused by the radial oscillation. The radial component of spin given 
by Eq. ()20j) is defined by the g— 2 precession. The importance of O has been first found by 
Granger and Ford and formula (|2~2~|) currently used has been derived by Farley 3|. 
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VI. DISCUSSION 



The particle EDM causes the spin rotation about the radial axis. As a result, the spin 
precession axis becomes tilted. If the spin rotation in the horizontal plane is cancelled 
or strongly restricted, the spin acquires a vertical component and the angle between the 
spin and the horizontal plane linearly increases from zero with time Q|. Therefore, the 
vertical magnetic field, the radial electric one, and magnetic focusing are expected to govern 
the particle motion in the EDM experiment. The best restriction of spin motion in the 
horizontal plane is obtained on condition that the radial electric field is adjusted to 

E = E = - — ^ B 

1 — ap 1 ^ 

sec Ref. 3). Such a restriction permits the high-precision detection of the EDM. 

Field misalignments and beam oscillations cause the spin motion imitating the presence 
of the EDM. In particular, the electric field may not be perfectly in the horizontal plane. 
In this case, the appearance of a vertical electric field leads to the vertical deflection of spin 
simulating the EDM effect . The use of Eq. (JT5|l simplifies the quantitative analysis of 
such effects. In this equation, the last but one term describes the spin rotation about the 
radial axis. The preceding terms also affect this rotation. The second-order effects discussed 
in Ref. [3J] can also be calculated by means of Eq. 1)1 5|) . This equation determines the 
general approximate relation between the effects of AMM and EDM: 

\n EDM \ _ _ P\d\ 
\u a \ 2\a\ \fj,'\ 

The EDM mig ht cause a small tilt of the spin precession axis of muons in the g— 2 



experiment 



18j . Measuring this tilt sets a limit of 2.8 x 10 19 e-cm on the muon EDM 



18| . If the muon EDM were equal to this limit, it would stimulate the rotation about the 
radial axis with the angular frequency IJI^da/I = 3 x 10 _3 |cj a |. 

The first excellent example of analytical description of complicated spin dynamics in 
storage rings is given in Refs. [f], Q, [a]. In these works, the effects of the vertical and 
radial CBOs on the horizontal motion of spin in the g— 2 experiment are described. For this 
purpose, the particle rest frame was used. However, this frame is inconvenient for taking into 
account the synchrotron motion and the particle acceleration caused by a locally nonzero 
longitudinal electric field. For example, such a particle motion is important in the EDM 
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experiment. If a strong radial electric field is used Q|, its imperfection results in the 
appearance of a longitudinal component leading to a systematical error. In Refs. j?l 1^ and 
other works, the effect described by the term oe z in Eq. (|14p was ignored. The description 
of spin motion given in these works is therefore approximate, while Eq. (I14j) is exact. 

Obtained equations (|1U| ) . (|14|) . and (J 15)) describing the spin dynamics in cylindrical co- 
ordinates have been derived in the general form. They can be helpful for high-precision 
experiments at storage rings when the spin motion needs to be analytically calculated. Nev- 
ertheless, these equations are convenient only when the ring is either circular or divided 
into circular sectors by empty spaces. Moreover, the analytical description given by these 
equations can be useless if the configuration of the main fields is not simple enough. In 
this case, it is necessary to use computer calculations. The algorithms of such calculations 
can handle all kinds of common magnet types and also the effects of misalignments. These 
algorithms are based on different coordinate systems (e.g., on Frenet-Serret coordinates \^). 
The method for calculating the spin motion in storage rings was developed by Derbenev, 
Kondratenko, and Skrinsky [l^| (see also Ref. and references therein). Another method 
using the Frenet-Serret curvilinear coordinates was proposed by Courant and Ruth These 
methods have been successfully used in various experiments and are also applicable in any 
other case. 

However, the analytical description of spin dynamics can be necessary for several high- 
precision experiments. Only the analytical description can guarantee the clear understanding 
needed for the g— 2 and EDM experiments. The use of Eqs. ([TUJ) . ([T4"Jl . and (fTK|) for analytical 
calculations of spin dynamics in these experiments can be very successful. Other formalisms 
seem to be less convenient. It is natural to use the cylindrical coordinate system for analyzing 
the spin motion in circular storage rings. Therefore, the obtained exact equation of spin 
motion in cylindrical coordinates can also be utilized in other cases where the configuration 
of the main fields is simple enough. 



VII. SUMMARY 



The exact formula for the angular velocity of particle motion in the horizontal plane has 
been obtained. The exact equation defining the spin motion in the cylindrical coordinate 
system with allowance for the particle EDM has been derived. This equation is convenient for 
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analytical calculations of spin dynamics when the configuration of the main fields is simple 
enough. Such calculations can be needed for several high-precision experiments. Cylindrical 
coordinates can be used if the ring is either circular or divided into circular sectors. Averaging 
of the angular velocity of spin rotation has been performed. The generalized formula for 
the influence of the vertical betatron oscillation on the angular velocity of spin rotation 
in the g— 2 experiment (pitch correction) has been found. This formula agrees with the 
previous result [?], 0] and contains the additional oscillatory term which can be used for 
fitting. The relative importance of the terms in the equation of spin motion including the 
EDM-dependent ones is discussed. 
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